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ABSTRACT 


A finite element analysis of hyperthermia in 
bioloqical tissues has been performed* The coupled thermal 
and electromagnetic phenomena have been studied in detailo 
The bio-heat transfer equation has been used for thermal 
modelling. Tv/o separate cases have been analysed, one 
with a fixed rate of blood perfusion through the tissue and 
the other with blood perfusion varying as a function of 
tissue average temperature and time. A two-dimensional 
problem with rectangular geometry has been considered for 
the sake of simplicity, although the formulation can be 
easily extended to 3-D or more complex geometries. The 
electromagnetic radiation is taken to be coherent and 
two plane polarized. 

The governing equation has been solved by the 
■finite element method using 8-noded isoparametric elements. 
Transient temperature results within the tissue have been 
obtained using the Frontal algorithm for matrix solution 
and marching in time using implicit finite differencing. 
Results have been presented for steady state temperature 
profiles, transient temperature variation at tissue centre 
and isotherm patterns, -^or various hyperthermia conditions. 
The effects of varying Specific Absorption Rate (SAR), rate 
of blood perfusion and the frequency of radiation, have been 

highlighted for two types of tissues, namely, the human 
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thigh muscle and the human lung. 

The numerical results have been validated by 
comparing with an analytical solution for the steady 
state oroblem vdth two adiabatic boundaries. The results 
indicate that the temperature increases vdth SAR linearly 
and decreases v;ith blood perfusion non-linearly. However, 
the temperature rise is less sensitive to changes in rate 
of blood perfusion. For the case of automatically changing 
blood perfusion rate, the tissue temperature displays an 
overshoot above the f^nal steady state value. At low 
radiation frequency, the rate of heating is more and it 
penetrates deeper. As for the electromagnetic properties of 
the tissues, the dielectric constant plays an important role 
in determining the level of temperature rise . 
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CHAPTER 1 
INTRODUCTION 

1 .1 GENERAL BACKGROUND ; 

Now-a-days cancer treatment is attracting more and 
more attention from medical personnels and technologists 
all over the world , as cancer is still one of the deadliest 
diseases which has defied man* Hyperthermia by high 
intensity microwave heating is one of the most effective 
methods that are being applied for cancer-treatment* The 
present work attempts to develop a detailed theoretical 
model for micro-wave heating of biological tissues, taking 
into account the electro*magnetic interaction between the 
tissues and the micro-wave radiation. 

It is well known that there is a critical temper- 
ature for biological cells (around 43°C - 45°C), beyond 
which the cell dies. This limiting temperature is called 
the tumoricidal temperature. If a certain portion of the 
human body containing the cancerous tissue is treated beyond 
the tumoricidal temperature, while the surrounding tissue 
is below this critical temperature, the cancerous cells can 
be destroyed selectively without affecting the other ceils 
of the body* 
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1 .2 DESCRIPTION OF HYPERTHERMIA ? 

In this process, the portion of the body, to 
be heated locally, is exposed to a high frequency electro- 
magnetic radiation. The frequency of the radiation may be 
in the short-wave radio-range or in the micro-wave range. 
Short-wave diathermy has been assigned the frequencies 
of 13.5, 27*1 2and 40.68 MHz, while the micro-wave diathermy 
has been assigned the frequencies of 915, 2450, 5850 
and 18000 MHz. ’Short-wave' diathermy provides deeper 
and more uniform heating than the diathermy at higher 
frequencies. Here we shall briefly discuss the inter- 
action between the electromagnetic radiation and the 
biological materials during hyperthermia. 

1^/hen the electro-magnetic wave passes through a 
certain medium, both electric and magnetic fields, which 
are perpendicular to each other and also to the direction 
of propagation of the wave, are produced. The interaction 
of these two fields gives rise to electromagnetic energy , 
which is carried by the travelling wave. When the wave 
passes through vacuum or air, the energy associated 
with it remains unattenuated with distance. But when 
the electromagnetic radiation passes through a dielectric 
medium of finite-electrical resistivity (such as the 
biological tissue), the energy of the EM radiation is 
attenuated gradually along the direction of propagation 
of the wave. This is due to the fact that the energy 
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of the wave is absorbed within the tissue at each point 
along the path of the wave. The electromagnetic energy 
absorbed by the tissue manifests itself in the form of 
heat, which raises the temperature of the tissue. The 
rate of heat absorption depends on the nature of the 
dielectric medium and also on the nature of the electro- 
magnetic radiation being incident* The rate of energy 
absorption increases for medium with higher values of 
dielectric constant. 

1 «3 DISCUSSION OF AVAILABLE LITERATURE ? 

Several electrical, medical and heat-transfer 
scientists have contributed to research on the analysis 
and technology of hyperthermia. 

Baish, Foster and Ayyaswamy [1 ] developed a 
parametric study for the relationship among the tube size, 
tube spacing, material properties and the simulated perfusion 
rate using perfused phantom models of microwave Irradiated 
tissues. They used a parallel tube heat exchanger config- 
uration to simulate the internal convection effects of blood 
flow. The measured thermal response of the ptiantom was 
seen to be favourably comparable with the numerical solution 
of the bio-heat transfer eguation under the same irradiation 
conditions. Thus they provided experimental evidence for 
the successful use of the bio-heat equation for modelling 
the thermal response of real tissues. 
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The same authors [2] developed analytical expre- 
ssions for the local temperature fluctuation near isolated 
and counter-current blood vessels during hyperthermia. 

These scaling lav;s relate the magnitude of such fluctuations 
to the size of the heated region and to the thermal equili- 
bration length of the vessels. They have shown that counter- 
current vessels have shorter equilibration length and produce 
smaller temperature fluctuations than isolated vessels of the 
Same size. Baish et al . [3 ] also conducted a parametric 
comparison of three different vascular mo dels_, namely, 

(i) an array of unidirectional vessels, (ii) an array of 
counter-current vessels and (iii) a set of large vessels 
feeding small vessels -which then drain into large vessels 
for describing heat-transport in tissue. Analytical and 
numerical methods were used to predict the gross temperature 
distribution throughout the tissue and the small scale 
temperature gradients associated with thermally significant 
blood vessels. They have sho-wn that three continuum formu- 
lations of bio-heat transfer (directed perfusion, effective 
conductivity, and a temperature dependent heat-sink) are 
limiting cases of the vascular models with respect to 
the thermal equilibration length of the vessels* 

V/einbaum, Jiji and Lemons [4] used the tissue 
clearance method for preparing vascular casts of -the rabbit 
thigh and sectioned them parallel to the skin surface to 
determine the detailed variation of the vascular geometry 
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as a function of tissue depth and proposed a nev three 
layer model for the deep tissue, skeletal muscle and 
cutaneous layers. In a latter work [5], they developed 
the above model into a detailed quantitative model which 
takes into consideration the variation of the number 
density, size and flow velocity of the counter -current 
arterio-venous vessels as a function of the depth from 
the skin surface, the directionality of blood perfusion 
in the transverse vessel layer and the superficial shunting 
of blood to the cutaneous layer. 

Weinbaum and Jiji [6 ]also developed a simplified 
3-D bio-heat equation to describe the effect of blood flow 
on blood-tissue heat-transfer. They derived a simple 
expression for the tensor conductivity of the tissue as a 
function of the local vascular geometry and the flow velocity 
in the thermally significant counter-current vessels. They 
have also shown that directed (as opposed to isotropic) 
blood-perfusion between the counter-current vessels can 
have a significant influence on heat-transfer in regions 
where the counter-current vessels are under 70- diameter. 

Arkin, Holmes, Chen [7] have performed a sensitivity 
analysis using a theoretical model of pulse decay method 
for the determination of local tissue thermal conductivity 
and blood perfusion rate. It is based on a comparison of 
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the measured temperatures with the theoretically calculated 
temperatures. They have shovm that a chosen pulse duration 
of 3 seconds is in aqreanent with the analysis as a good 
compromise between accuracy and excessive tissue heating. 

The same authors have described in a later paper 8 , 

the theoretical model upon which the TPD method is based and 

have detailed its capabilities and limitations. 

Valvano, Allen and Bowman [9] have presented an 

improved technique for the determination of the thermal 

conductivity, the thermal diffusivity, and the blood 

perfusion using a self -heated spherical thermistor probe. 

In this paper, it has been shown that in the presence of 

—1 /2 

flow, the transient power response does not follow t 
as has been previously assumed. Here, firstly, effective 
thermal properties are measured and then perfusion is quan- 
tified from effective thermal properties. 

Sekins, Emery, Lehmann and |»iacDougall [10] have 
shed light on the blood flow response occurring during 
local hyperthermia in muscle through a two-dimensional 
transient thermal model of human thighs undergoing 
microwave diathermy. The model blood flow values were 
checked against those measured in the human subjects via 
xenon 133 wash out method and good agreement was found-. 

Foster et al . [11 ] have examined the heat- 

transfer characteristics of biological bodies which are 
subjected to strong microwave fields with simultaneous 
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cooling by flowing water to estimate the maximum temperature 
increase and the thermal tissue constants that are encountered 
durin g typical hyperthermia experiments. 

Oilier and Hayes [12] have modelled the burn process 
resulting from the aoplication of a hot, cylindrical source 
at the skin surface using the finite element technique. 

Here a rotationally symmetric 125 element mesh was defined 
within the tissue beneath and outside of an applied heating 
disk. Natural convection with ambient air was assumed for 
areas of skin surface not in direct contact with the disk. 

Nev/man and Lele [13 ] have proposed several techniques 
to measure the local thermal diffusion in tissue by heating 
with a focussed ultrasound field. Transient as well as 
near steady state heat inputs have been considered and 
examined for their suitability of being the measurement techni- 
ques for either tissue thermal diffusivity or perfusion 
rate. 

Elkowitz, Shitzer and Eberhart [ 14 ] have numerically 
solved the bio-heat transfer equation to calculate the tern- 
perabure profiles in tissues subjected to non-uniform blood 
flow distribution for initial and boundary conditions which 
simulate experimental physiological situations. The methods 
and results are useful for the prediction of temperature 
profiles in the absence of significant endogenous or 
exogenous heating. 
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J.F, Deford and O.P.Gandhi [15 ] have used the impedance 

method for solving low frequency dosimetry problans in 
three-dimensional geometries involving lossy dielectrics* 

The scattering object has been modelled with an equivalent 
impedence network. Their results indicate that the contact 
electrodes produce relatively uniform deposition between 
the plates in the region of torso, whereas the loop appli- 
cator studied gives rather large deposition near the surface 
o-^ the body, especially in the spinal region, and very small 
deposition in the core of the body* 

J. Behari [16] has dealt with the phenomenon of 
dielectric dispersion and mechanism of microwave interaction 
in biological media* He has also cited several electrical 
methods for measuring the electrical properties of tissue* 

O.P. Gandhi, J.Y* Chen and A. Rlazi [ 17 ] investi- 
gated how currents are induced in a human being subjected to 
microv/ave radiation. They have also discussed about the 
safety qu3.delines regarding the induced current* 

Thus the available literature Indicates that the 
research studies have so far concentrated selectively on 
the physiological or thermal or electrical aspects of 
hyperthermia in biological materials* More research is 
required to be performed on the coupled nature of the 
various processes during hyperthermia* 
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■> -4 SCOPE AMD OBJECTIVES OF THE PRESENT STUDY : 

The present study deals with a coupled thermal 
and electromagnetic model of hyperthermia in biological 
tissues by microwave radiation. For the sake of simplicity, 
a two dimensional rectangular geometry has been considered. 
However, the problem formulation is of a very general nature 
and it can easily be extended to axi-symmetric and three- 
dimensional problems. The applied electromagnetic radiation 
is assumed to be polarised such that the wave propagates only 
in x-direction. 

The objectives of the present work are to predict 
the temperature distribution within a tissue of the human- 
body and ascertain the portions which will be destroyed by 
radiation. The following specific aspects of the problem 
will be considered! 

i) To predict the range of power of the micro- 
wave radiation that can be used to destroy 
the cancerous tissue selectively. 

ii ) To show the effects of different radiation 
frequencies upon the temperature profile at 
the same power input and blood perfusion 
rate for a given tissue. 

iii) To demonstrate the effects of constant and 
time dependent blood perfusion rates upon 
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the heat-transfer characteristics of the 
hyperthermia process. For the time 
dependent case, the blood perfusion rate 
vail be assumed to depend upon the average 
tissue temperature and also explicitly upon 
time . 

iv) To highlight the influence of the electro- 

magnetic properties of the tissue by performing 
the analysis for different tissue materials. 
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CHAPTER 2 

thermal modelling 


2.1 INTRODUCTION : 

In the earlier days, the detailed processes which 
occur during the passage of microwave radiation through 
biological tissues were not well understood. It was 
believed that the destruction of tissues by microwave radi- 
ationwas due to electrical ohenomena. Later it became 
established that thermal effects chiefly contribute to the 
cell destruction during hyperthermia. In the present study 
of hyperthermia by microwave radiation, the electro-magnetic 
interactions are talten into account through Maxwell's field 
equations and the thermal processes are modelled via the 
bio-heat equation. 

Many studies have so far been performed on heat- 
transfer analysis in biological tissues using the bio-heat 
equation. But none of than have so far used the Maxwell's 
field equations to predict the rate of heat generation for 
unit volume of tissue. In these analyses, a uniform heat 
source of given intensity has been assumed. In the present 
study, an attempt has been made to predict the rate of 
heat generation at any given location using only the 
electro-magnetic properties of the tissue material. and the 
characteristics of the microwave radiation. Thus actual 
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problem parameters such as the intensity of radiation 
and the microwave frequency are directly involved in 
the analysis. Also the specific nature of the tissue 
being subjected to hyperthermia is taken into account 
through the use of appropriate thermal properties (i.e* 
density, specific heat and thermal conductivity) and the 
electro-magnetic properties such as the electrical resis- 
tivity, electrical permittivity and magnetic permeability. 

Another important parameter in bio— heat equation, 
which controls the rate of tissue heat-up, is the local 
blood perfusion rate. The proper evaluation of this 
parameter is complicated by the fact that the physical 
and mental resoonse of the human body to any stimulus, 
say^that of the micro-wave heating, is very complex. It 
has been observed that the blood perfusion rate is different 
at different locations of the human body. It also varies 
with time, temperature of the tissue, physical and 
mental conditions of the subject and so many other environ- 
mental factors. An accurate correlation of blood perfusion 
rate with its influencing factors is yet to be established. 
Keeping all these facts in mind the present study is 
restricted to two specific situations of blood perfusion. 

i) In the first case, the blood perfusion rate 
has been considered to be constant with time and tissue 


temperature. 
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ii) In the second case, blood perfusion rate has 
been considered to vary with time and the average temper- 
ature of the tissue. An experimental correlation -From [ 1 O] 
has been used to calculate the blood perfusion rate. 


2.2 MATHEMATICAL MODELLING : 

The transient heating of a rectangular tissue with 
microwave radiation has been considered. A two-dimensional [Fig. 2.1 
situation has been modelled here for the sake of simplicity 
although the formulation can be very easily extended to 3-D 
or axisymmetric geometries*. The electromagnetic radiation 
is taken to be corresponding to that of a polarized uni- 
directional wave. The rate of blood perfusion is assumed 
to be uniform throughout the tissue. The initial temperature 
of the tissue at time t = 0 (just at the start of heating) 
is taken to be equal to the normal body temperature of 
36°C. The outside atmospheric temperature is assumed to 
be 27°C. The thermal and electromagnetic properties [Table 2.1] oi 
the tissue are taken to be uniform spatially and constant 
with respect to time. With these assumptions, the temper- 
ature-time response of the tissue has been estimated using 
the bio-heat equationo 

2.2a BIO-HEAT EQUATION ; 

The energy conservation principle applied to the 
biological tissue undergoing diathermic or microwave heating. 


Convection 



1cm 

J_ ::: 

EM 

Radiation 


FIG. 2.1 TISSUE REGION C( 



T= 36 C 



6 cm 


ED FOR SOLUTION DOMAIN 


Table 2.1 ! Value of the Thermal, Electrical and 

Magnetic Properties for Air, Muscle, 
Lung and Blood 


For Air or Vacuum ^ 


•"o = 

4 X 10"’^ 

9 9 

N-S /C 

Yo = 

0.8854x0“^^ 

C^/M-m^ 

i^or Muscle 

and Lunq: 


^Pt 

= 2.925.0 

J/kg-°C 

P 

t 

= 1000,0 

Kg/m^ 

'<t 

= 0.625 

W/m-°C 


= 1 .0 


r 



For Muscle 

only? 


^ r 

= 128 


r 

= 1 .5873 Ohm-ffi 


e 



For Lung only: 


Yr 

= 40 



r = 3.7037 Ohm-m 

e 

T^or Blood t 

3 

P , = 1 000 .0 kg/m 

b 

C = 4184.0 J/kg-°C 
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may be expressed by the transient bio-heat transfer equation 
of the forms 

K. v^T - W. C , (T - r ) + Q = P. C , (2.1 ) 

^ bpb a tptat 

where T is the local tissue temperature, is the tissue 
thermal conductivity, Wj^ is the blood perfusion rate, 

Cpb is the specific heat of blood, is the arterial blood 
temperature (generally considered as constant), Q is 
the rate of volumetric heat-generation, Cp^ is the effective 
specific heat of the tissue, is the density of tissue 

and t is the time of heating. 

The above formulation ^or bio-heat transfer 
analysis is chosen, because it qualitatively describes the 
tissue characteristics reasonably well and it provides an 
adjustable parameter, Wj^, that has been identified with 
the rate of blood perfusion in a real tissue. The blood 
perfusion term in the above equation acts as a heat sink, 
which arises from thermal equilibration of the blood in 
the capillaries with the surrounding tissue. Recent studies, 
that rigorously relate the blood tissue heat-transfer to 
the vascular organisation in tissues, have shown that blood 
thermally equilibrates in relatively large counter-current 
vessels [5,6] * They suggest that the bio-heat equation 
seems to work only because it provides enough adjustable 
parameters to fit the available data. In any case, for want 
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of a better model at present to describe tissue hyperthermia, 
the bio-heat equation has been adopted in the present work. 

2 ,2b INTERACTION OF EM RADIATION WITH DIELECTRIC MEDIA : 

The Maxwell’s equations which govern the variation 
of the magnetic a*nd electric field vectors (H and E respect- 
ively) are given by: 

ae 

2 X H = Y * 3^ 

a H 

V X E = - p • 

V . E = 0 (2.4) 

V , H = 0 (2o5) 

The properties of electromagnetic waves through any dielectric 
medium are governed by the Maxwell’s equations given above. 

For simplicity, the following restrictions are imposed on 
the electromagnetic radiation to which the biological tissue 

i s expo sed : 

(i)) A plane wave of incident radiation is consi- 
dered that is propagating in x-direction. All the properties 
associated with the wave are assumed to be constant over 
any YZ plane at any given instant of time. Thusj 


( 2 . 2 ) 

(2.3) 


ay 


az 


0 



18 


(ii) The electro-magnetic wave is polarized such 
that the electric field vector E is contained only in the 
X-Y plane and the magnetic field vector H is contained only 
in the X-Z plane. Thus, 

= 0 , E^ = 0 

Hy = 0 , = 0 


After incorporating the above-mentioned assumptions 
in equations (2.2) - (2.5), the general electric and magnetic 
field equations are obtained ast 




and 



W • 











( 2 . 6 ) 


(2.7) 


Equations (2.6) and (2.7) can be identified as the wave 
equations that describe the propagation of the electrical 
and magnetic component in the x-direction for an isotropic 
medium of finite electrical conductivity. 


Analytical solution for equation (2.6) iss 


= Eyj^ [ exp i6)(t “ x)] exp (- ~ KX) 


( 2 . 8 )/ 


and that for equation (2.7) iss 


H- 


exp [i<o(t - X)] exp (• 


Q 


KX ) (2.9) 
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By observing the above two equations, it is quite obvious 

that Ey and are attenuated exponentially with distance 

when the radiation passes through an isotropic medium of 

finite conductivity. The term exp(- KX) in both 

o 

the equations (2.8) and (2.9) causes the attenuation of 
the electromagnetic field intensities. The complex refra- 
ctive index, n, for the dielectric medium is 'given bys 

n = + n - iK 


where K is termed the extinction coefficient for that 
medium as it causes attenuation of the electro-magnetic 
wave. The values of n and K can be found out from 
the relations: 






[1 + [1 + (■ 




1/2 


( 2 . 1 0 ) 


2 itC, 


r Y 
e‘ 




W c 


[-1 + [ 1 + (■ 


2 1/2 

•) ] 3 


2tc C 


r Y 
e 


( 2 . 11 ) 


It can be shown that the relation between Ey and is 

(n - iK) 




uc. 


Ey 


( 2 . 12 ) 


2.2c ENERGY FLUX DETERMINATION OF JHE EM WAVE : 


The instantaneous energy carried per unit time 
and per unit area by an electromagnetic wave is given by 
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the cross-product of the electric and the magnetic field 
intensity vectors. This product is known as the Poynting 
vector S, given by the expression! 

S = E X H (2.13) 

and according to the properties of the cross-product, S 
is a vector propagating at right angles to the electric 
and the magnetic field vectors in a direction defined by 
the right hand rule. For the present case, the magnitude 
of S is given as 


I SI = direction of S is the 

same as that of the wave propagation. 

Using equation (2.12) in (2.13), 

ISI = Ey (2.14) 

0 

Thus the instantaneous energy per unit time per unit area 
carried by the vrnve is proportional to the square of 
the amplitude of electric field intensity. 

Considering the real part of Ey only. 




(2.15) 


Substituting for Ey equation (2.14) from equation (2.15), 


n 


pCo ^YM 


exp(- 


SwKX 


) f(X,t) 


(2.16a) 
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where, 

f(X,t) 


j^Cos [oo ( t - 


n 2 

X)]] 


(2.16b) 


For the present case, since the electromagnetic 
radiation is of very high frequency, the instantaneous energy 
can be averaged over time ast 


2w KX 


* ^ ^av. 


wh ere. 


f (X, t) 
av . 


2 Vo) 

2i7<0 f 


Thus, 


S I = 


n 


E 2 u KX 

^ exp[- [ ] ] 

2 Cq 


{2.17a) 


(2.1 7b> 


( 2.1 8 ) 


Considering the real part of equation (2.18), 


n 

pC 


2m 


'o 


2 


exp [- 


2wKX 


] 


(2.19) 


The extinction coefficient exists only for an 
electromagnetic wave passing through an isotropic medium 
of finite electrical conductivity. For the electromagnetic 
wave passing through air or vacuum, since these media can 
be assumed to possess infinite electrical resistivity 
fi.e. r = on), the extinction coefficient K is zero and 
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honce there is no attenuation term in the analytical solution 
for the electric and the magnetic field intensities. In 
such a case: 


^Y,0 [w(t - ^ X) ] 

and 

0 ~ ^ZM (t - -Q X) ] 

f o 


(2*20a) 


(2.20b) 


where E and ^ are the field intensities in air ©r 
Y, 0 4, u 

vacuum. The instantaneous flux of an electro— magnetic x^rave ^ 


while passing through air, is then given by 




(CONSTANT) 


( 2.21 ) 


For the present problem, ISJ is nothing but the 

..mO 

energy of the incident electromagnetic radiation before 
entering the tissue region of interest. 

2 .2d CALCULATION OF THE RATE OF HEAT GENERATION PER UNIT 
VOLlffAE: 

As the electromagnetic wave penetrates the 
biological tissue, due to the inherent dielectric proper- 
ties of the tissue medium, some portion of the energy 
carried by the wave is absorbed within the tissue. This 
rate of energy absorption per unit volume of the tissue 
is given by: 
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Q = - (v . s) 


(2o22) 


Since S exists only in x-direction , 

3 3 2oKX 

Q = - ^ (ISO = - exp(-— ) ] (2.23) 

The negative sign in the above expression arises due to 
the fact that | S l decreases in the positive x-direction* 

Thus, 

Q = 4 • »P <- ^2.24) 

11 Co ™ % 

Using equation (2*21) in equation (2*24), 



His 

U 


exp 




(2.25) 


vjhere, 

Q = 
o 



It is noted from equation (2.25) that the volumetric 

rate of heat-generation depends upon so many factors like 

Q , n,K, \i etc. Q is a dimensional parameter indepen- 
0 * r o 

dent of the di-electric medium and is dependent only upon 


the energy and the frequency of the incoming wave. 
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2 *26 GOVERNING EQUATIONS ? 

After substitution of Q from equation (2*25) the 
bio~heat equation (2*1) can be revnritten as follows; 


2 

V T 


«b Sb - V' + % <- 


2(0KX X 
C ^ 
o 



aj[ 

at 


( 2 . 26 ) 


Boundary Conditions t 

Here, the heat transfer analysis is performed over 
a region v/hich is sufficiently large compared to the 
dimension of the Irradiated area. The boundary conditions 
that have been used for the four faces of the solution 
domain C^ig. 2.1) are those of convective heat loss to the 
atmosphere and known boundary temperatures* For the face 
which is exposed to the atmosphere, convective boundary 
condition has been applied in the forms 

Kt lx ^ h(T - T^)at X=0and0<Y<L (2.27a) 

Since the portion , which experiences the effect of hyper- 

thermia, is small compared to the whole region, the three 
faces of the region at Y = 0, Y = L and X = L can be 
considered to be at the body temperature. Thus, 
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X=L, 0 < Y L for all t 

T=T^for Y=0, 0 ^ X < L for all t 

Y=L, 0 < X < L for all t 

Initial Condition ! 

Prior to the exposure to microwave heating, the 
entire tissue region is assumed to be at normal arterial 
temperature. Thus, 

T = T3 

at t=:0 for 0 < X < L and 0 < Y ^ L (2«27e) 

2.3 DETERMIMATION OF BLOOD PERFUSION RATE ? 

In the present study, the blood perfusion rate 
has been estimated in the following manner: 

(i) Firstly, the blood perfusion rate has been 

considered to be constant with time and unaffected by 

the tissue temperature. In this case, blood perfusion rate 

3 3 

has been varied from 1 kg/m s to 7<.0 kg/m s, independent 
of the SAR (Specific Absorption Rate) value. 

I 

(ii) In the second approach, the blood perfusion 
rate has been considered to vary with time and tissue 
temperature. This assumption is much more realistic in 
nature. The temperature control systan inside the human 
body has a very complex behaviour. The body tries to restore 


(2.27b) 
(2. 27c) 
(2.27d) 
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the temperature to the normal value of 36°c, whenever 
it decreases or increases. This is done partly by 
varying the blood perfusion rate and partly by other 
mechanisms such as perspiration. Generally blood perfusion 
rate increases with the increase in tissue temperature. 

This control of blood perfusion rate for different tissue 
temperature is done by the brain through biological feed- 
back control mechanism. There is always some time-lag 
for the signal to reach the proper destination. Hence, 
blood perfusion rate is also a function of time. 

The nature of dependence of blood perfusion rate 
on time and temperature has been taken from [10]- Here 
simple linear interpolation and extrapolation have been 
used to find out the blood perfusion rate, corresponding to 
a given average tissue temperature at a given time from 
the data presented in Table 2.2. These blood perfusion 
data For human thigh muscles of persons having an average 
build, are shown in Table 2«2o 


2.4 NON-DIME NSIONALISATION OF THE GOVERNING EQUATI.ONj.t 


The bio-heat equation has 


using the following dimensionless 


2 ? 

- , I 


Y 

L ’ 


been non-dimensionalised 
variables s 

" - H 

— T — T 

a amb 


* 

t 




Table 2.2 i Variation of Blood Perfusion Rate v/ith Tissue Temperature and Time for 
Human Thigh Muscle [10 J 
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Equation (2.26), can then be rewritten ass 


V T + Q exp(-BX ) - W T = 

a t 


(2.28) 


where. 


Q 




^t ^'^a”'^anib ^ 


nK 


(2.28a) 


W 


vr 

b „ _ 2<^KL 


(2.28b) 


and, 


2 0) S, 


(2.28c) 


The non-dimensional initial and boundary conditions ares 


ii} 

(ii) 


T =0 at t =0 for 0 < ^ < 1 and 0 < Y < 1 (2.29a) 

3 T* hL ^ ^ "^f “ "’'a ^ r * * 

^ [T - Tf - i t : :) J ° - Tfi 


3X 


arab ' 


at X =0 and 0 <Y < 1 


(2«29b) 


-9f # ^ 

(iii) T =0 at X =1 and 0 < Y < 1 


(iv) T =0 at Y =0 and 0 < X < 1 


(2.29c) 

(2.29d)/ 


(v) T =0 at Y =1 and 0 1 X < 1 


(2.29e) 
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CHAPTER 3 

finite element analysis 

3.1 INTRODUCTION: 

The solution procedure that has been used for 
the present study is the finite element method. There are 
many advantages of using FEM as the solution procedure for 
the hyperthermia analysis. First of all, it can handle any 
complex' geometrical shape of the tissue being irradiated. 

The heating or cooling boundary conditions for the tissue 
of interest may be introduced in a very general manner. 

Material property variation with location can be very 
easily accounted for. In view of the above mentioned 
generality and flexibility provided by FEM solution technique , 
it has been employed for solving the present problem. 

F *2 DESCRIPTION OF THE FEM SOLUTION PROCEDURE BY GALERKIN 
APPROACH : 

In FEM, the solution domain vhich is continuous, 
is divided into several subregions each of which is known 
as an element. Each of these elements is of chosen shape 
and a prescribed number of nodes are placed within and on 
the boundary o-f each element. The collection of elements 
within the entire solution domain is called the finite [Fig.3.1 3 
element mesh. The solution variables are approximated 
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within each element by suitable approximating functions. 

'^^e error in satisfying the governing equations is minimised 
in an integral sense over the whole domain leading to the 
formulation of a set of simultaneous algebric equations for 
all the nodal variables. Finally, the solution of these 
nodal equations gives the required solution. The above steps 
can be summarised as follows! 

(a) The equivalent integral equations are obtained 
from the governing differential equations by minimising the 
error in satisfying the governing equations in an integral 
sense over the whole domain. 

(b) The continuous solution domain is divided 
into a number of elements v/ith a chosen number of nodes 
placed at desired locations within each element. All the 
elements can be of the same shape or of different shapes. 

A mesh-generation program may be used to obtain the finite- 
elementHmesh [Fig.3.l]. 

(c) ) Each governing integral equation over the 
solution domain is written as a sum of the integral over 
all the elements. 

(d) A typical element is isolated and the approxi- 
mating functions for describing the field variables are obtained 
within that element in terms of their nodal values. For 

this element, firstly, the elemental integrals ar® evaluated 
leading to the Formation of an elemental coefficient matrix 
and elemental right hand side vector. Then the boundary 



32 


integrals (if any); are evaluated to obtain the boundary 
coefficient matrix and the boundary right hand side 
vector* These elemental and boundary coefficient matrices 
assembled for each element to obtain the final form of 
the elemental coe-f^f icient matrix* A similar operation is 
done for the elemental right hand side vector also. 

(e) The elemental coeTficient matrices are then 
assembled for all the elements, giving rise to a global 
coefficient matrix. The elemental right hand side vectors 
are also assembled leading to the formation of the global 
right hand side vector. 

(f) The matrix equations obtained in the previous 
step are solved to find the unknown field variables at all 
the nodes. 

3.3 METHOD OF WEIGHTED RESIDUALS ; 

Weighted residual methods are, in essence, 
numerical techniques which can be used to solve a set 
of partial differertial equations. The first step in 
the application of the weighted residual procedure is to 
assume that the unknown variable, 0, can be approximated 
over the whole domain by 
* 

0 = S N, 0. 

1=1 ^ ^ 


(3.1) 
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where are the interpolation functions described in 
terms of the spatial coordinates (X, Y) and 0^ a»« the 
values of 0 at certain chosen locationso 

Consider the equation: 

^((0> = 0 (3.2) 

vtiere 0 is the exact solution of the differential operator 
An approximate solution, 0 , given by eg. (3*1) gives rise 
to an error or residue R while substituted into eq. (3*2), 
such that 



In the weighted residual methods R is minimised 
over the whole solution domain in the integral sense with 
a set of arbitrary 'weighting functions', such that 

over the whole domain, D, 

/ W- R d D = 0 (3.3) 

D ^ 

where i = 1,2, 

If the number of unknown values of 0 is n, then 
by using n linear independent weighting functions 
enough equations can be generated to solve all the unknowns. 
The only limitation, at this stage, plaped on is that 

they must be positive, single-valued and finite. 
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The Galerkin’s Method is a very special case of 
the method of weighted residuals, where the weighting 
functions, V/^, are taken to be the same as the interpolation 
functions, N^, which are more commonly known as shape 
functions . 


Therefore, for Galerkin 's Method, the residue 
minimization equations can be written as: 



The dimensionless bio-heat-transfer equation (with 
the stars omitted for convenience) is given as* 

v^T - WT + Q exp (-BX) = (3.5) 

Applying the Galerkin method to this equation, 

we obtain the residue minimisation equation as: 

a T 2 "BX 

If ^ dXdY -If V T dX dY - // R Q e dX dY 
A A A 

+ // N, WT dX dY = 

A 


0 


(3,6) 


35 


Thp unknovm variable, T, can be expanded in 
terms of its nodal values and shape functions as; 

T = S N. T. (3,7) 

j=1 ^ ^ 

where n is the number of nodes per element. 

Using equation (3.7) in equation (3.6) and applying 
integration by parts to the second term of eqn. (3.6), 
we obtain* 

[srt\ Nj dx dY] tTj 1 + [ /n^ (-k^ dl ] 

t 3 

3N. 3N. 3N. 3N,. 

+ ^ dX dY] } 

- [Q // dX dY] + [W // dX dY ] { Tj 1 

= 0 (3»8) 

^“-^t boundary contribution 

to the whole domain. 

In the above equation, the summation, S , has been 
dropped by adopting the convention that summation is implied 
if a subscript is repeated. 

In the present study, 8-noded iso-parametric 
elements have been used [See Fig. 3.1] . For isoparametric 
elements, the global coordinates (X, Y) of a point inside 

I 

t' 

t, 
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an element can be expressed in terms of nodal coordinates 

and shape functions as it has been done for the field 
variables. Therefore, 

8 

^ = 2 NX (3.9^. 

i=1 

8 

^ = 2 NY (3^.1 0> 

i_1 -L i 

The shape functions Nj_ for the 8-noded iso-parametric 
element in terms of local coordinates | and t) , are 
given by the following relations. 

For corner nodes: 

\ = ^(l+£j^£)(l+h Tl)(^.g + Tl T1_1) 

1 1 i 

(3«11) 

For mid-nodes: 

Nj_ = for =0 (3*12) 

N. = “ (1 -Ti^){l +£ §) for h. = 0 (3.13) 

1 ^ 1 1 

Substituting equations (3*9) - (3o13) in 
equation (3*8);, the final form of the assembled matrix 


equations are: 
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[GMASS] ir. ] + [GC0N]{TJ = {GGEN} (3,14) 

where [ G^'ASS] is the global transient matrix, [ GCON ] is 
the global coefficient matrix,! GGEN} is the global right 
hand side vector* 

In terms of shape functions or their derivatives, 
these matrices and vector are expressed as! 


m 


[GMASS] = 2 N.n. dX dY 

e=1 (e) ^ J 


(3.15) 


[GCON] 


m 


S [W // 
e=1 Ce) 


a N. aN. a N, aw. 

N^Nj dX dY + (— . -g-a + 

dXd' 


+ C2 ^ N, N. dl ] 
(e) ^ J 


(3.16) 


m 


(GGEN ) = E [ Q e' 

e=1 (e) 


'BX 


N^ dXdY + C1 N^dl ] 


with Cl 


■q for heat flux condition 


*■ 


and C2 


Bi for convective B.C. 


0 for heat flux condition 


(3.17) 
(3 .18a) 

(3.18b) 

(3ol9a) 


Bi for convective B.C. 


(3.19b) 


where q* 
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Here the line integrals v;ith coefficients C1 and C2 arise 
due to boundary conditions. 

The equation (3.8) represents a set of simultaneous 
first-order ordinary differential equations of the variable 
T and its time derivatives. 

Here the finite difference technique has been adopted 
to approximate the time derivatives# This leads to the final 
matrix equations of the forms 


[A] iTJ = {C} 


(3 .20) 


There are three major schemes for handling the time 
derivatives in the finite-difference technique, namely: 

i) Implicit scheme 

ii) Explicit scheme 

iii) Crank -Nich olson scheme. 

Applying a generalised scheme -which considers a 
fraction of the explicit scheme end e fraction of the implicit 
scheme, eq . (3.14) can be rewritten as, 

K+1 

[[GMASS] + 0^ t [GCON]] t Tj! 


= [[GMASS] - (1 - 0)a t [ GCON]]tTj1^ 

+ { GGENlAt 


(3.21 ) 


The three schemes mentioned above can be obtained 
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from eqn. (3.21) by assuming different values of the 
■fraction. For instance, 

^ = 0 for explicit scheme 

0 = 0.5 for Crank-Nicolson scheme 

' © =1.0 for implicit scheme. 

Since the implicit scheme is the most stable one 
and has no restriction as far as the time step is concerned* 
it has been used in the present study. 

’=’or implicit scheme, equation (3.211 simplifies to 

the form: 

[[GMASS ] + A t[ GCON]] 

K 

= [ GMASS] t Tj I + { GGEN ! A t (3 .22 ) 

Equation (3.22) can be solved to obtain the temper- 
ature vector at the (K+1)th time level, knowing the temper- 
ature vector at the Kth time level. 

3.5 MATRIX SOLUTION PROCEDURE S 

Here the Frontal solution method has been used to 
solve the assembled matrix equations. It is a well-known 
fact that the method adopted for solving the assembled 
matrix equation has a significant bearing on the computer 
storage requirement and execution time when the total number 



40 


of unknov/n variables is very large. It becomes very time 
consuming and the core memory problem is very accute v/hen 
a large number of unknowns are solved by the matrix 
inversion technique. The Frontal solution technique can 
handle quite a large number of variables without the problem 
of large core memory and the computer execution time is 
also qui te less. 

The Frontal solution technique used in the present 
X s 

work/based on the direct Gaussian elimination procedure for 
solving the symmetric matrices where the leading diagonal 
is always used as a pivot. For unsymmetric matrices encoun- 
tered in a wide range of engineering problemSt the most 
suitable oivot is not necessarily on the leading diagonal 
and a Frontal solution routine exists for off diagonal pivoting 
[ 20] . But since this tends to be more time consuming, 
the method used in the present work uses only diagonal 
pivoting and incorporates many features of the Frontal 
method for solving symmetric matrices. 

The overall Frontal solution technique covers 
the following steps t 

^ Poj^ation of the elemental matrices, 

il) Assembling Into a small, temporary global matrix 

of chosen size. 
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(ili) Introduction of known variable boundary condition* 

(iv) Reduction of the qlobal matrix usinq Gaussian 
elimination procedure. 

(v) Back substitution. 

The primary objective of the Frontal method is 
the elimination of variables as soon as possible after their 
introduction, via appropriate equations, into the qlobal matrix. 
When the contributions from all the elements to a particular 
node point have been assembled, the corresponding variables 
associated with that node can be eliminated. The complete 
matrix is , therefore, never assembled, since the reduced equations 
can be eliminated from the core and stored on disc* The 
equations held in core, with the corresponding nodes and vari- 
ables, are termed as the FRONT and the number of unknown 
variables in the front is termed as the FRONT WIDTH. 

The Front width changes continuously, since, if 
the nodal matrix contribution have been fully summed the 
corresponding equation reduction based on the diagonal 
pivot can be executed. 

For a non-symmetrio global matrix, a preassigned 
global matrix core ar« la filled from contributing elements, 
the largest diagonal entry in the pre-assigned core is then 
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found and used as the pivot in a direct Gaussian elimination 
process o As the maximum, predetermined number of equations 
93:e eliminated, the corresponding reduced equations are 
written on to disc and more elements and corresponding 
equations are taken into core. 

The equations, nodes and variables currently in 
core are termed active, those assigned to disc as deactivated 
and those yet to appear in core as inactive. This is shown 
diagrammatically in Fig. (3.2). 




O Inactiv® NodiS 
• Active Nod«s 
O Deactivated Nodes 


mZA Front 


Next Element of Assembly 
t t : Assembled Element 


Fig. 3.^ Basic idea of Frontal method. 
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CHAPTER 4 

RE.SULTS AND DISCUSSIONS 


Transient temperature results have been obtained 
for hyperthermia of human thigh muscle and lung tissue 
for various SAR values, blood perfusion rates and frequency 
of radiation. The results are presented in the form of 
steady state temperature variations with depth, transient 
temperature variation at the tissue centre and isotherm 
patterns v/ithin the tissue. For checking the validity of 
numerical results, the special case of hyperthermia with 
two adiabatic boundaries (Fig. 4.1 ) and one convective and 
one prescribed temperature boundaries havebeen considered. 
Analytical solution exists for the steady state temperature 
distribution for this case (Appendix A). 

In ^ig. 4.2, the steady state temperature variation 
with depth along the centre-line is compared for the 
analytical and numerical solutions at SAR = 1.0x10 W/m , 
=5.0 kq/ms and «= 27 MHz. The t^^ results agree 
very clearly (within 0.059^), indicating the validity of 
the present numerical result by FEM# 

In Figs. 4.3a, 4.3b and 4.3c, the steady state 
temperature profiles within the muscle tissue are shown 
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var.ious combinations of SAR and V/, , Due to the symmetry 

D 

of the boundary conditions used in the present formulation, 
the results are identical about the mid plane of Y = 0 . 

For this reason, in all the three figures, only teiTiperature 
variation of one half the tissue is showoo Initially , the 
temperature increases slightly with depth and then decreases 
gradually and finally drops rapidly to the prescribed boundary 
temperature of 36'^C. The reason for the maxintum temperature 
not occurring at the surface but slightly into the interior 
is that there is some loss of heat to the atmosphere due to 
convection at the outer boundary which is exposed to the 
atmosphere. The gradual decline in temperature with depth 
is partly because of the decrease in the intensity of the 
EM wave and partly because of the heat conduction in the 
lateral direction. The sudden drop in temperature near 
the inner boundary indicates that the radiation is not 
fully absorbed due to the finiteness of the domain considered. 
A much larger depth is required to completely dissipate all 
the EM energy for the parameters assumed for these figures. 

The figures also indicate that the maximum temperature is 
achieved along the centre line of the applied radiation, 
as expected. There is a sharp decrease in temperature 
in the lateral direction away from the mid-plane and for 
distances roughly the order of radiation beam diameter, 
the temperature^As°aimLt%lose to the normal body temperature 
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Similar trends are observed for all combinations of SAR 
and V/[^ , except that for high blood perfusion or low SAR 
value, there is a small region far away from the tissue 
mid-plane, along the exposed face where the temperature 
IS below 36 C. This is because of the convective cooling 
on the outer boundary. 

In Figure 4.4, the temperature profiles along 
the tissue mid-plane are presented for small time and 
large time for the thigh muscle. Effects of various 
combinations of SAR and Wj^ have also been plotted. It 
is seen that, at all times, the temperature profiles are 
similar to the steady state profiles discussed earlier. 

The graph also indicate that for very large value of SAR 
or small value of blood perfusion, the temperature levels 
are high and vice-versa. However, the temperature change 
is more sensitive to SAR variation than it is to 
variation. The time taken to reach steady state is found 
to be higher at larger SAR or lower Wj^'fhis due to the 
fact that a significant portion of the absorbed energy is 
carried away by the blood when the steady state is reached. 
At lower blood perfusion rates, due to the lower heat 
carrying capacity of the blood, it takes longer time to 
reach steady state. At higher SAR also, more amount of 
heat has to be removed by the blood , which takes longer 


time. 
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The transient variation of muscle tissue centre 
temperature is showi in fig, 4*5 for various SAR values 
at fixed The steady state value of temperature 

increases linearly with SAR. This feature is corroborated 
by the analytical solution of appendix A, The time taken 
to achieve steady state slightly increases with SAR as 
discussed earlier. 

In fig. 4.6, the variation of muscle tissue centre 
temperature is plotted with time for different Wj^ values. 

The temperature decreases with an increase in blood perfusion 
as heat sink effect is higher at larger Wj^ . The steady 

state temperature (at the tissue centre) varies approximately 

- 0.6 

as Wj^ . The analytical solution also shows that 

the dependence of steady state temperatures in the middle 
portion of the tissue will fall between l/VWj^ and 1 • 
Because of this, temperature variation is not as sensitive 
to variation in blood perfusion as it is to the variation 
in SAR value. 

The effect of the EM frequency upon the temperature 
rise and the heat generated is shown in figs. 4.7a and 4. 7b. 
The surface temperature is essentially insensitive to the 
radiation frequency, whereas the interior temperature is very 
much influenced by the frequency of radiation at lower 
frequency levels. The temperature profile is essentially 
seen to be governed by the rate of volumetric heat generation 
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within the tissue (Piqo 4«7b), except near the inner and 
outer surfaces ,virh ere the applied boundary conditions play 
an important role. Since the rate of volumetric heat 
generation decays very steeply with depth and lower amount 
of heat is generated at high frequencies, the temperature 
values are also lower and display a strong decayo 

The comparison between the mid-plane temperature 
profiles has been provided for the muscle tissue and the 
lung tissue in fig* 4*8* The figure shows that the temper- 
ature profiles are qualitatively similar for both types of 
tissue. For the lung tissue, the temperature levels are 
lower chiefly because of lower dielectric constant. All 
the other important electric and magnetic properties such 
as permeability, resistivity etc. are of the same order 
for both the tissues. It is seen that the decay constant 
is also lower for the lung tissue. 

In Fig. 4.9, the transient tenperature variation 
at the tissue centre is shown for the muscle and the lung 
tissue. The increase of temperature is seen to be 
roughly proportional to the dielectric constant. The larger 
the dielectric constant, the more is the dissipation of energy. 
The dependence of heat generation on dielectric constant is 
seen to be linear. The time taken to reach the steady state 
is larger for the muscle. 
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figure 4.10 shows the temperature profiles along 
the tissue mid-plane for small time and large time for 
the lung. It essentially shows all the features mentioned 
for '^ig« 4.4 except for the range of temperature variation. 
The range of temperature variation for lung tissue is much 
less. This is due to the low value of dielectric constant 
for lung. 

The temperature variation i^rith time at muscle tissue 
centre is plotted in Figo 4.11 for automatically varying 
blood perfusion rate at two different SAP. levels. It is 
seen that the temperature increases rapidly in the beginning 
and reaches a maximum value at around 18 minutes after the 
start of heating and then decreases with time. The initial 
rapid increase in temperature is because of low rates of 
blood perfusion. However, as the tissue heats up, the 
temperature-control mechanism is initiated by the brain 
and the rate of blood perfusion starts increasing with time. 
As a consequence, the rate of increase of temperature within 
the tissue is slowed down. After a typical response time 
of the order of 15 minutes, the temperature starts decreasing 
due to the already increased value of blood perfusion. 

The temperature then has a tendency to settle dovm at a 
much lower steady state value. It is also seen from the 
graph that with higher SAP value, the initial rate of 
temperature increase and the later rate of temperature 
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decrease are both steeper. The steeper drop of temperature 
beyond the maximum value occurs due to the larger value 
of at high temperaturee 

In Fig* 4.12, the temperature-time response at 
the muscle tissue centre has been shown for varying blood 
perfusion and constant blood perfusion models, at a fixed 
SAR. The fixed model shows a monotonic temperature 
variation tnv/ards the steady state value. This is because 
of th® fact that the heat-sink effect arising from blood 
perfusion increases linearly v/ith tissue temperature, 
and at steady state, an equilibrium situation is reached 
when a significant portion of the supplied heat is removed 
by blood perfusion. The varying model, on the other 
hand, has non-linear heat-sink behaviour due to the 
dependence of on tissue temperature. Since there is 

also time lag between the rate of temperature increase and the 
rate of Wj^ increase, the temperature shoots above the final steady 
state value. 

In Fig. 4.13, the muscle temperature profile for 
varying are presented for different SAR values. For 
the sake of comparison, temperature profile with constant 
W,^ is also shown. It is evident from the graphs that 
the temperature profile is similar for fixed or varying 


W. at all SAR levels. 
D 
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In figures 4.14a,b,c,d, the isotherm plots within 
the muscle tissue are shown for fixed and varying at 
two difFerent time levels. The isotherms in the heated 
egion are U-shaped loops indicating that the heat absorbed 
is mainly conducted away in the lateral direction. The 
temperature gradient in the axial direction is not very 
high near the exposed surface. For larger depths, however, 
the effects of the applied boundary conditions comes into 
play so that the isotherms near the inner boundary are 
closely spaced. This features also indicates that all of 
the applied radiation has not been absorbed within the 
prescribed depth in the present problem. Far away from the 
heated region, at some locations near the surface exposed 
to atmosphere ^ temperature falls below 36°C because of 
the convective heat loss to the atmosphere. This consequence 
occurs because the metabolic heat generation has been 
neglected in the present work. A comparison of figures 
4.14a and 4.14b shows that the temperature profiles are 
similar at large and small time levels except that heat 
penetrates to a larger distance for larger irime. In the 
case of varying at small times, heat penetration is 
confined to a small region. With increasing time, the ■ 
entire tissue is heated and some portion of the tissue gets 
overheated as the variation of blood perfusion rate is still 
lagging behind. At a still larger time, sufficient increase 
in blood perfusion has taken place so that the tissue is 
again cooled. This phenomenon is supported by the 

Fig » ( 4 .1 4e) e 



b2 

Thus the numerical results provide both qualitative 
and quanhitative picture of the effects of various para- 
meters durinq hyperthermia by micro wave radiationo 
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FIG. 4-5 TRANSIENT VARIATION OF MUSCLE CENTRE TEMPERATURE 
AT VARIOUS POWER LEVELS 
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FIG, 4 9 A COMPARISON OF TRANSIENT TISSUE CENTRE TEMPERA- 
TURE FOR MUSCLE AND LUNG 
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FIG. 4. 11 MUSCLE CENTRE TEMPERATURE VARIATION WITH TIME FOR 

VARYING BLOOD PERFUSION RATE 
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FIG. 4.Ua ISOTHERM PATTERN IN MUSCLE FOR FIXED 
BLOOD PERFUSSION RATE 
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FIG. 4 .UC ISOTHERM PATTERN IN MUSCLE FOR VARYING 
BLOOD PERFUSION RATE 
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FIG.4.M ISOTHERM PATTERN IN MUSCLE FOR VARYING 
BLOOD PERFUSION RATE 


FIG 4.Ue ISOTHERM PATTERN IN MUSCLE FOR 
VARYING BLOOD PERFUSION RATE 
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CHAPTER 5 

CONCLUSIONS AND SUGGESTIONS 

5#1 Conclusions 

(1) Detailed numerical results have been obtained 
whose validity has been checked with an available analytical 
solution for a special case. 

( 2 ) The increase in temperature level is more 
sensitive to 3AR than to blood perfusion rate. 

( 3 ) Vlhen blood perfusion rate varies automatically 
(as in a real tissue) the temperature has an overshoot and 
subsequent decrease to the final steady state value. This 
is due to the time lag associated with the increase in blood 
perfusion rate, as a response to the Increase in tissue 

temperature. 

( 4 ) At lower frequencies, the volumetric rate of 
heat generation is higher compared to that at higher frequen- 
cies. Also, the decay of the heat generated with depth, 

is less steeo for lower frequency levels. Therefore at low 
frequencies, the temperature level is higher and it decays 

more slowly with depth. 

5 .2 Suggestions 

( 1 ) The present problem can easily be extended to 
three dimensional or other complex geometries. 



(2) The phenomenon of variable blood perfusion 
rate can be look^ in ho more deeply* 

( 3 ) More complex situations, which considers th 
effect of (i) unidirectional and countercurrent blood 
vessels, (ii) micrnstructure of tissue, (iii) metabolic 
heat p'^neration etc*, can be focussed deeply for the 
present studyo 
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APPENDIX A 


An analytical solution of the bio-heat transfer 
equation (2.28) is possible for steady state ''dth boundary 
conditions as shov/n in Pig, ( 4 .I). For such a situation, 
the problem becomes one-dimensional and the non-dimensional 
■joverninq differential equation is: 



VIT + Q = 0 


(A.1 ) 


To obtain the solution of equation (A.1), let us assume 


T = Tk(^)+ T (X) 
h p 


(A. 2 ) 


where Tp(X) is the particular solution and Tj^(X) is the 
homogeneous solution, we shall assume Tp{X) to be of 
the form: 


T 


P 



^-BX » "3X 

e = C e 


(A. 3) 


where, 



Substituting for T in equation (A.1>, we get. 


17 “ 


W = 


0 


(Ao4) 
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The general solution of equation (A, 4) is given 

by : 

T}^ =: Sinh VW X + Cosh/W X (A. 5 ) 

where and are constants to be determined from boundary 
conditions. The complete solution to equation (A. 1 ) is, 
there'f^ n re * 

T = A-] SinhVW ^ + A^ CoshsTV/ X 4 C e"®^ (A. 6 > 



This boundary condition gives rise to: 

A = y [ ^2 + + -Ii' - A 3 

V Vi 

ii) At X = 1 , T = 0 (A. 8 a) 

This boundary condition gives rise to: 

A^ SinhVW + A^ Cosh/ W + C e"® = 0 (A. 8 b) 


Solving equation (A. 7 b) and equation (A, 8 b), A^ 


and A^ 


can be determined as: ^ „ 

-Ce-®+ [y -CO + |r)] SlnlVW 



/TT 


SinhV W + Cosh / W 


(a *93) 
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A 


1 


/ V/ 



* 


^ C(1 + 



(A. 9b) 


Substitutim ^or 
equation (A. 6) gives the 


the constants A^ and A^ in 
complete solution for temperature. 



